Quantum criticality analysis by finite size scaling and exponential basis sets 
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We combine the finite size scaling method with the meshfree spectral method to calculate quantum 
critical parameters for a given Hamiltonian. The basic idea is to expand the exact wave function 
in a finite exponential basis set and extrapolate the information about system criticality from a 
finite basis to the infinite basis set limit. The used exponential basis set -though chosen intuitively- 
allows handling a very wide range of exponential decay rates and calculating multiple eigenvalues 
simultaneously. As a benchmark system to illustrate the combined approach, we choose the Hulthen 
potential. The results show that the method is very accurate and converges faster when compared 
with other basis functions. The approach is general and can be extended to examine near threshold 
phenomena for atomic and molecular systems based on even-tempered exponential and Gaussian 
basis functions. 



I. INTRODUCTION 

The study of how the energy levels of a given system 
change as one varies a parameter in the corresponding 
Hamiltonian is of general interest, particularly near bind- 
ing threshold, level crossings and quantum phase transi- 
tions. In phase transitions, critical points are associated 
with singularities of the free energy which occur only in 
the thermodynamic limit [ll, Q- Finite size scaling (FSS) 
was developed by Fisher and others ^^ to calculate 
such parameters by extrapolating information from a fi- 
nite system to the thermodynamic limit. In analogy, FSS 
was also developed to extrapolate information from a fi- 
nite basis set to the infinite basis set limit in order to 
calculate quantum critical parameters for a given Hamil- 
tonian. This is done by expanding the exact wave func- 
tion in a complete basis set and use the number of basis 
function to play the role of system size Q. Early work 
using FSS to calculate quantum critical parameters was 
based on expanding the wave function in Slater-type and 
Gaussian- type functions p3l-[llj. Recently, the method 
was also combined with the finite element method (FEM) 
13, D3 and B-splines expansion to achieve similar results 




Here, we combine FSS method with the meshfree spec- 
tral method (SM) to calculate quantum critical param- 
eters. Lately, the meshfree SMs start gaining growing 
attention because of their high levels of analyticity and 
accuracy |15l - [l9| . In these methods, the unknown func- 
tions are approximated by expansion using preselected 
basis sets. One of the main challenges in SM is to handle 
domains extended to infinity [20. - 26J . Many techniques 
were introduced to overcome this challenge such as us- 
ing exponentially decaying functions as basis sets, the 
truncation of the computational windows, and applying 
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size scaling. Recently, a non-orthogonal predefined ex- 
ponential basis set for eignevalue problems involving half 
bounded domains was introduced and used [27, '28'| . The 
set is easy to use and allows generally finding a wide 
range of eigenvalues simultaneously. 

In this paper, the exponential basis sets are imple- 
mented in FSS analysis to obtain the quantum critical 
parameters for a given Hamiltonian. The presented tech- 
nique is real-space meshfree. Such real-space techniques 
start gaining more attentions in "ab initio" and density 
functional calculations [l2, l20| . As a benchmark system 
we choose the Hulthen potential. For such Hamiltonian, 
the analytical solution is known and also FSS was imple- 
mented using other basis functions, hence our numerical 
results can be compared and analyzed. The comparison 
confirms the validity and efficiency of the new approach 
and its applicability for FSS analysis which will be used 
on more complex systems. 



II. THEORETICAL BACKGROUND 

A. Analytical solution for the Hamiltonian with 
Hulthen potential 

Hultheii potential |30l . |31| is a special case of Eckart 
potential [32| , which is a family of screened Coulomb po- 
tentials. It has the following form 
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where A is the coupling constant and a is the scaling 
parameter. For small r, it resembles Coulomb potential. 
But, it dies faster and exponentially for large r. 

By defining a dimensionless variable, x — -, and in- 
serting the potential in Schrodinger radial differential 
equation, the radial equation becomes 
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The analytical solution for this equation is known in term 
of hypergeometric function and it is 

^ = iVoe""^' (l - e-^) 

2F1 (2a + l + n, l-n,2a+l;e-^) (3) 

where a = —a^E, n is the state order, and iVo is the 
normalization factor and it is 



iVo = ^/a [a + n) (2a + n) 



The energy levels are 
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It is clear from Eq.(IS]) that Ac = n? /2 is a critical 
coupling constant. As Ac is a function of n, it is obvious 
that the number of allowed bounded states (i.e. Umax) is 
A dependent. It has at least one state for A > 1/2 and 
this is the critical point to be tracked. 



B. Finite size scaling 



where v is the scaling exponent for the correlation length. 
To obtain the numerical values of the critical parameters 
(Ac, a) for the energy, we define for any given operator 
O the function 
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If we take the operator O to be i? ~ Eth and dH/dX, 
we can obtain the critical parameters from the following 
function [8] 
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which at the critical point is independent of N and N' 
and takes the value of a. Namely, for A = Ac and any 
values of N and N' we have 



r„(Ac,iV,7V') = a- 



(11) 



Because our results are asymptotic for large values of 
N, we obtain a sequence of pseudocritical parameters 
(Aat, a^v) that converge to (Ac, a) for TV — > cx). 



As aforementioned, FSS method is a systematic ap- 
proach allowing extrapolating the critical behaviour of 
an infinite system by analysing a finite sample of it. It 
is efficient and accurate for the calculation of critical pa- 
rameters of the Schrdinger equation. Assuming that the 
Hamiltonian of a system is of the following form 



H = Ho + Vx (A) 



(6) 



where again A is the coupling constant. The critical 
point, Ac, will be defined as a point for which a bound 
state becomes absorbed or degenerate with a continuum. 
As known, the asymptotic behaviours of physical quan- 
tities near to the critical points are associated with crit- 
ical exponents. So, the energy near to Ac can be defined 
as 
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where we assume that the threshold energy, Eth does not 
depend on A. In principle, Ac can be calculated providing 
the exact solution. However, when use variational calcu- 
lations to expand the exact wave function of the system 
in a basis set, only a finite number of basis functions [N] 
can be used practically. So, the calculated physical ob- 
servable (i.e. Ex in this case) depends on N . Thus, for 

each N^ the calculated energy level is denoted by _Ej^ ' . 
FSS assumes the existence of a scaling function Fe such 
that 
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C. Spectral methods and the exponential basis sets 

Meshfree SM is a special family of the weighted resid- 
ual methods |15l - [l§ |. In these methods, the unknown 
functions are approximated by either an expansion of or 
interpolation (known as collocation method) using pres- 
elected basis sets. For homogeneous and smooth compu- 
tational windows, SMs work very well. But, they suffer 
from the Gibbs phenomenon if any of the structural func- 
tions of the studied problem is not analytical. To avoid 
this problem, the computational window is divided into 
homogeneous domains where the discontinuities lie at the 
boundaries. This approach is known as multi domain 
spectral method (MDSM) [il-[l|. In general MDSM 
methods allows handling very complicated and discontin- 
uous functions. This capability is very flexible as any ex- 
pansion basis set can be used. In this paper, the studied 
problem has a smooth structural function (i.e. Hulthen 
potential). So, MDSM is not used. 

In many physical problems, the extensions toward in- 
finities decay exponentially as 
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where ± is used to cover both =Fcx) with positive /3. As 
aforementioned, this is one of the main challenges in SM 
[20l - [26j . A review paper by Shen and Wang discusses 
this problem in further details [2a|. Recently, a non- 
orthogonal predefined exponential basis set for eignevalue 
problems involving half bounded domains was reintro- 
duced [23, Sal- Similar sets were introduced in 1970s 



by Raffenetti, Bardo, and Ruedenberg |33l - l35j for self- where /?„ is as defined above in Eq. [Ml 
consistent field wavefunctions. 

The set is easy to use and it overcomes many chal- 
lenges such as zero-crossing and single scaling problems 
by approximating the decaying domain functions by ex- 
ponential basis set which spans wide range of decaying 
rates as follows: 



III. IMPLEMENTATION 



Formulation 
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where c„ are the expansion coefficients and /3„ are the 
pre-selected decaying rates. They are chosen intuitively 
based on the studied problem. But, they should allow 
many possible decay rates with very small number of 
bases. In this paper, the decaying rates are defined as 



/3„ = 10^" 
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where ds and de are the smallest and largest used powers 
respectively and N is the number of the used bases. 

In this paper, the set is modified slightly to have a 
faster convergence by enforcing the states to vanish at 
X = 0. The modified set is 



To simplify the moments calculations, the normalized 
Schrodinger radial differential equation (Eq. [2]) is rewrit- 
ten as 
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The expansion form (Eq. [T6l) is used to solve the above 
equation. For each used number of basis iV, the expan- 
sion form is rewritten as follow: 
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This form is working only for bounded states and hence 
should work fine only for A > 0.5. By implementing this 
expansion form, Eq. [T7]can be written as: 
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where the elements of the matrices are the following 
scalar products: 
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and by selecting proper values for dg and de , a wide range 
of eigenvalues {E\^ J ) and their corresponding eigenstates 

(c3^ ) can be calculated directly. The used values for dg 
and de are -4 and 4 respectively. 
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In the above three equations, the integrations are taking 
place in one dimension and not over the physical three 
dimensional space. Eq. [TO] is a direct eigenvalue problem 



In this paper, we focus on the critical change in the 
lowest energy level. So in the remaining of this paper, 
Ej^ J is corresponding to the calculated ground state level 
with N basis; clearly, it is a function for A and the scaling 
factor a. Also, c\^ ' is corresponding to the ground state 
and it contains the expansion coefficients. Generally, the 
states need normalization by dividing the coefficients by 
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In this case and the foUowing calculations for the po- 
tential energy, the integrations are calculated over the 
physical three dimensional space for the case of Z = 0. 

To apply FSS as shown later, we need to calculate the 
potential energy. It is simply 



yr = -4.A^cW*cW 
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(24) 
The integrations are computed numerically by Gaussian 
quadrature. Obviously, this is the most numerically ex- 
pensive part in work. However, it is clear also that the 
integrations are independent of the state distinctive pa- 
rameters (i.e. A and Cn )■ So, for each N^ the integra- 
tions are calculated at the beginning and the results are 
used to calculate V^^ while varying A. 

To obtain the critical parameters, we use the following 
shifted functions: 
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The critical parameters Ac and a can be obtained from 
the Fq, as defined in Eq. [TOl 



B. Results and discussion 

In the calculations, the scaling parameter (a) is set to 
one. Also as aforementioned, the used parameters for the 
exponential basis set are -4 and 4 for ds and dg respec- 
tively. These parameters are chosen after few iterations 
to have a reasonable accuracy for the eigenvalues. To im- 
plement FSS, N is varied between 32 and 48 in a step of 
2. So, Ac and a can be obtained by seeking the crossing 
of the FSS curves. 

The calculated ground state energies {Eq) are shown in 
Fig. [T]as a function of A for all the used values of N. The 
errors are very small (as shown in Fig. ^ and hence the 
lines are overlapping. More resolution (in A) is shown in 



the small box. As can be observed, the calculated values 
for the ground state energy start diverging slightly from 
the exact solutions as A approaches Ac- This is expected 
as the used basis works for bounded states and the error 
shall increase as the states get extended in space. How- 
ever, the calculated values of Eq are still very accurate 
and a relative error of about 10~^" was obtained around 
A = 0.51 for iV = 32 and A = 0.5001 for iV = 48 as shown 
in Fig. H 
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FIG. 1. The calculated ground state energy (-Bo) as a func- 
tion of A using different numbers of bases, which are varied 
between 32 to 48 in steps of 2. The errors are very small (as 
shown in Fig. [2]) and hence the lines are overlapping and thus 
more resolution about Ac is shown in the small box. 
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FIG. 2. The relative errors of the calculated ground state 
energy as a function of A for A*' = 32 and A'^ = 48. 

In Fig. 131 the results of FSS calculations are shown. 
Plotting Fa as a function of A for different values of N 
gives a family of curves that intersect around the ana- 
lytical Ac — 0.5 and a = 2. The exact crossing of any 
adjacent curves defines the pseudo-critical parameters A^r 
and a AT, which are used to analyse the convergence. 
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FIG. 3. Fq as a function of A. The numbers of bases are 
varied between 32 to 48 in steps of 2. 
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FIG. 4. The extrapolated values of Ac vs. 1/A'^ as obtained 
by this work (solid circles), Hermite interpolation polynomials 
(HIP) (white squares [i3])j and finite element method (FEM) 
(crosses Q2 



To check the convergence, the pseudo-critical parame- 
ters Xn (Fig. HI) and un (Fig. [5]) are plotted as functions 
of 1/N and compared with the results obtained using 
Hermite inter pola tion polynomials (HIP) FEM by An- 
tillon et al. [l3|. It is clear that the three methods 
converges to the analytical values. However, the used 
exponential basis set in this paper results in consider- 
ably faster convergence when compared with the other 
two methods. The results of the three methods are sum- 
marized in Table I. 

The last point to be presented is to confirm the va- 
lidity of FSS assumptions using data collapse calcula- 
tion. In Fig. ini EqN""/" is plotted as a function of 
(A — Ac) N^^/" for all the used N values. It is clear that 
all the curves overlap perfectly and thus validates our 
FSS assumptions. 



FIG. 5. The extrapolated values of a vs. 1/N as obtained by 
this work (solid circles), Hermite interpolation polynomials 
(HIP) (white squares [13|1), and finite element method (FEM) 
(crosses [13l[). 



TABLE I. Results for critical parameters. 
Analytical This work FEM fl3] HIP fl3] 



Ac 
a 



0.5 
2 
1 



0.500001 
2.00094 
1.00000 



0.50184 
1.99993 
1.00079 



0.50000 
2.00011 
1.00032 



IV. CONCLUSION 

In atomic and molecular physics, the near threshold 
binding is important in the study of ionization of atoms 
and molecules, molecular dissociation and scattering col- 
lisions. Our benchmark calculations for the near thresh- 
old behavior of the energy levels of the Hultthen po- 
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FIG. 6. Data collapse study of the used approach using dif- 
ferent numbers of bases, which are varied between 32 to 48 in 
steps of 2. 



tential indicate the validity of combining FSS method 
with the meshfree SMs to calculate quantum critical pa- 
rameters. Fortunately, the exponential basis sets used in 
this study have been used previously as exponential-type 
even-tempered basis for atomic orbitals :33-35i] . The re- 
sults indicate that cvcn-tcmpered bases are very accurate 
in Hartree-Fock atomic calculations. Also, a system- 
atic approach extending even-tempered atomic orbitals 
to optimal even-tempered Gaussian primitives have been 
developed and used decades ago in standard quantum 
chemistry calculations for atomic and molecular system 
[si IM IM, ill. Thus, our combined FSS method and 



SMs based on even-tempered basis sets might be used 
to extract quantum critical parameters for atomic and 
molecular systems. In future studies, we plan to combine 
our FSS procedure with the Hartree-Fock and density 
functional theory (DFT) and other "ab initio" methods 
using SMs with even-tempered basis and other intuitive 
basis sets to analyse criticality and near threshold phe- 
nomena for molecular and extended systems. The pre- 
sented approach allows scaling to analyze large systems 
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